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Abstract. In this paper an effective quantum particle beam-splitter 
in the momentum space is realized in the frame of open-loop control 
scheme. We demonstrate for small interaction time that the splitting ef- 
fect ±40ftfc with summarized relative intensity in both main components 
is about 50 per cent from initial intensity of the atomic beam. 
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1. Introduction 

In nano-lithography the optical control of atomic motion is one of the 
main problems. In principle, we can make such a control for atom dynam- 
ics because there is an exchange of momentum between the atoms and the 
optical fields. The momentum exchange can be created with practical de- 
vices: atomic mirrors and atomic beam splitters which are main elements 
of atomic interferometer. The most interesting possibility here is to ob- 
tain the splitting of the initial atomic wave packet coherently into two main 
momentum components only by controllable way. It is needed both for 
increasing of atomic interferometer sensitivity and for the creation of pe- 
riodic nano-structures by atomic wave packet lithography Previously, 
successful splitting of an atomic wave packet has been achieved by using Ra- 
man pulses, magneto-optical beam splitter, diffraction in an optical standing 
wave, adiabatic passage. More recently, coherent splitting has been realized 
by scattering of an atomic wave packet in standing wave with modulated 
intensity [2] and by using chirped standing wave fields 3j. 
In this paper we concentrate on the possibility to split an atomic wave packet 
in standing wave with modulated amplitude because this beam-splitter has 
a number of advantages by comparing with others. The first one is the sim- 
plicity for an experimental realization because it is quite easy to obtain the 
time modulation of intensity with any shape. The second advantage is that 
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the scale splitting of an atomic wave packet can be controlled by chang- 
ing the values of both an amplitude and the frequency of the modulation. 
Actually, we demonstrate for small interaction time, which is required for 
clear and large splitting, the beam-splitting effect zizAOhk with summarized 
relative intensity in both main components about 50 per cent from initial 
intensity of the atomic beam. 

The principal opportunity to split the beam in the momentum space was 
demonstrated in [3]. To achieve the effective splitting, we apply the scheme 
of open-loop control, or feed-forward control, i.e. a control signal depends 
only on the time. Our control goal is to obtain the large angle splitting for 
the initial wave packet after some time of the interaction between the atoms 
and the field of modulated standing wave. 

2. Physical background and mathematical model for beam 
splitting in the momentum space 

Atom lithography is an active field now a days. The resolution of an 
optical lithography technology is limited by diffraction, which for the case 
of deep ultraviolet light approaches 200 nm. The progress of recent device 
technology requires smaller patterning of 10 nm size. However, when one 
tries to make very small devices, the resolution of the resistance is limited 
by the spread of the secondary electron in an electron beam lithography, as 
well as in X-ray lithography. 

The ability to generate ultracold atoms using lasers has opened up new 
possibilities. The long de Broglie wave of cold atoms makes possible an 
interferometric manipulation with atomic wave packets, which is designed 
by an optical standing wave. In this case, atoms can be controlled directly 
to form a desired pattern. To produce the pattern with high resolution, 
we need to split the wave packet into two coherent momentum components 
only. For the model with only two states (i.e. an approximation of two level 
atom states), we have to split the population of the lower state (because 
the population of the exited state usually loses the coherency very fast by 
spontaneous decay) in several momentum components (in an ideal case - 
only two). At the same time to form the pattern with small step we need 
to control the scale of splitting between two main coherent components in 
momentum space. Therefore, an atomic beam splitter is the main element 
for the practical realization of nano-scale lithography with the controlled 
step by coherent scattering of an atomic wave packet. 

Let's consider now a two level atom in a far detuned standing wave with 
the intensity modulated in time as I = Iof{e,A ■ t,(j)Q)cos{kx), where e 
is an amplitude, A is the frequency of the modulation, and (pQ is an initial 
modulation phase. The standing wave with the frequency ivi applies between 
two states of atom system, where the state 1 is ground and the state 2 is the 
exited one. Here, loq is the frequency of atom transition and the difference 
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LOi — ujQ is the detuning. We will assume that the beam from an atom 
source propagates along z-axis and crosses the optical wave, standing along 
X-axis, by right angle. The spontaneous emission from the upper level in this 
system can not be neglected. After some time t of the interactions between 
the atoms and the field of the standing wave, the initial atomic wave packet 
is splitted in few coherent momentum components. 

Dynamics of the atom in the modulated standing wave is described with 
nonstationary Schrodinger equation for the wave function '^{r,t) of the two 
level atom: 

ih^^ = mr,t), (1) 

where is a Hamiltonian which takes into account both the atom movement 

along the standing wave and the dipolc interaction between the atom and 
the optical field. For sufficiently large detuning, when it is much larger than 
Rabi frequency and the natural width of the atomic transition, J7 >> Ro,T 
(where i?o is the Rabi frequency, T is the natural width of the transition) , the 
excited state 2 can be adiabatically eliminated. As the result, we obtain the 
equation for the amplitude of the probability of the ground state 

in ' ' = -— A,,*i(ar, t) + -^[f{e, A • t, ^o)?cos-'{kx)^i{x, t), (2) 

and m is the atom mass. 

After the Fourier transform the same equation in the momentum space is 
given by: 

= + R'')^iip, t) + ^ [^,ip + 2, r) + - 2, r)] , (3) 

where = ^^^^ [/(g; {^/^r) • t, 0o)]^ , = fik'^/2m is called a recoil 
frequency, r = LORt. Here we normalized atom momentum along x-axis to 
hk and the other quantities (the interaction time, the Rabi frequency and 
the detuning) we normalized to the recoil frequency lur. We have to point 
out that equations (2)- (3) are valid in the approximation when both the 
changing of atom momentum along z-axis and the initial value along x-axis 
can be neglected. 



3. Shell model for the splitting process: parametric control 

To explain the effect of splitting in the momentum space we start from 
the case of parametric control with a constant R. We invent a complex shell 
model for ^'i-function. Initially the atomic beam has a Gaussian distribution 
centered at p = 0. Thus, from the structure of the RHS (3) we can expect the 
non-zero meanings of functions to be concentrated in the neighborhoods 
of the points p = 2n, where n = 0, ±1,±2,.... Then we can predict the 
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continuous dependency on p with the discrete number n: + 2n,r) = 

y„(r). Dynamics Eq. (3) can be re- written in the form: 

= (4n2 + R')yn{T) + ^ [yn-iir) + Vn+iir)] (4) 

with the initial conditions: ?/o(0) = 1 y„-^o(0) = 0. Now we want to Umit 
our shell number. In the case of three shells only Eq. (4) becomes very 
simple: 

= R'yoir) + R'yiir); (5) 

i^ = {R' + ^)yiir) + ^yo{T). 

We demand for the cider shells : y±2 = y±4 = ... = for any moment r. Eq. 
(5) can be easily solved: 

yo(T) = e-^(^'+2)^ [Cie"^^ + C2e-"^^] ; (6) 
yi(T) = -le-^(«'+2)- [-{00 - 2)Cie'-- + (a; + 2)C2e-^--] , 




with uj = \/2R^ + 16/2. For the initial conditions yo(0) = 1 and ?/i(0) = 
the constants are: C\ = {io + 2)/2uj and C2 = {uj — 2)/2uj. The corresponding 
population amplitudes of the shells and ±1 are given by: 

, ^ , ^ *, ^ -, ■ 2I V2i?4 + 16r\ 
ao(T) = yo{T)yo{T) = 1 - „^ , ^ sm ( ) ; (7) 

ai(r) = yi(r)yt(r) - ^ 

Surely, the normalization ag + O-i + a-|_i =09 + 2ai = 1 is saved for any 
moment r. Now we can see that in the 3-shell model the regular splitting 

is realized when the time t = {2k + l)7r/\/2R^ + 16, fc = 0, 1, 2, The 

same effect can be reproduced in the case of five shells. To simplify the 
final expression we will omit the coefficients 4 and 16 in RHS, because the 
numerical meaning of is about 300 (i.e. i?^ >> 4 and 16). Then we 
apply Laplace transform: 

i{sYo{s)-l) = R^Yo{s) + R^Yiis) 

isYi{s) = R''Yi{s) + ^[Yo{s) + Y2{s)] (8) 

isY2{s) = R^Y2{s) + ^Yi{s). 
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Yois) 



Yiis) 



IE 



?2 



(9) 



Y2is) = - 



With the inverse Laplace transform we restore the time-dependent solution 
in the momentum space: 

iR^T \ f iR^T 



yo{T) = 



exp{—iR t) + exp I — 



2(2 + V3) 



+ exp - 



2(2 - V3) 



2/1 W = 



exp 



iR^T 



2(2 - \/3) 
y2{T) = -^exp(-ii?^T) + ^ exp 
The amplitudes are harmonical: 

ao{T) = yo{T)yo{T) = 



exp 



iRh 



2(2 + V3) 



(10) 



2(2 + V3) 



+ exp 



2(2 - V3) 



1 + 2cos y^R^r 



aiir) = yi(r)yt(r) = -sin^ —R't 



(11) 



a2('r) = y2{T)yl{T) 



l_cos l^^i^V 



Surely, again the normalization ao+a_i+a+i+a_2+o+2 = ao+2ai+2a2 = 1 
is satisfied. The splitting effect is obtained, when cos{\/?>R^t /2) = —1, then 
the population of the ±2 shells is: 2a2 = 8/9, and in the same time ai = 
and oo = 1/9 only (see Fig.l). 

However, if the number n of a shell is increased such that Am? » R? (i.e. 
for R ~ \/300 we have n» 10), then R^ in (4) can be excluded as a small 
parameter, and for the elder shells 

■dynir) 



dr 



c± An^yn{T) {n» 10) 



(12) 



This function is almost independent of the neighbor shells and it has the 
solution 

yn{T) ^ e-^*"'-y„(0) . (13) 
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Figure 1. Splitting effect for the model of 5 shells. 

But yn(0) = for any n ^ 0, thus, the elder shells do not participate in the 
re-distribution of the initial Gaussian population. Thus, the simple para- 
metric control with the fixed R is not enough to split the beam efficiently. 
Another scheme of time-dependent R (corresponding to the most general 
open-loop control) should be applied. 



4. Numerical simulation results for open-loop control with 

harmonical modulation 

Now let us consider the two level atom in a far detuned standing wave 
with an intensity, which is modulated in time harmonically as / = /o(l + 
ecos(At))^cos^(A;j;), where e is the amplitude and A is the frequency of the 
modulation. 

We assume also that an initial wave function r = 0) has Gaussian 
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profile with tlie widtli 5p: 



0) 



exp 



{5pf 



(14) 



We remind that now we use the dimensionless time r = uoRt. Fig. 2 shows the 
numerical solution of an equation for amplitude of the probability of ground 
state t)P in momentum representation for the cases unmodulated 

and modulated standing wave. We assume that initial wave packet has the 
width equals 6p = O.bhk and e = 0.8, A/w/j = 29. As we can see from this 
picture, the scattering result strongly depends on the amplitude modulation 
exiting in this system. If for an unmodulated case, it is well-known scatter- 
ing picture observed (Fig. 2), when an initial wave packet is splitted into a 
number of momentum components. However, for modulated standing wave 
the scattering picture is changing dramatically and two main momentum 
components centered on zizAOhk can be observed (Fig. 3). Such behavior of 
the momentum components is due to specific parametric resonance, which 
occurs in this system by the well defined amplitude and frequency modu- 
lation. We have to point out that the values of the modulation obtained 
for the amplitude and frequency modulation are strongly different from [2j , 
and we can interpret such resonances as Bragg resonances of high orders in 
modulated standing wave. 

I- 
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Figure 2. The dependence of the distribution function on 
an atom momentum for an interaction time Tint = 0.567. 
The unmodulated standing wave with the dimensionless Rabi 
frequency Rq = (320)^/2. 
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Figure 3. The dependence of the distribution function on 
an atom momentum for an interaction time Tint = 0.567. 
The modulated standing wave with the dimensionless Rabi 
frequency Rq = (280) 
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